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Abstract 

An explicit expression for the numbers A(n, r; 3) describing the refined 3-enumeration 
of alternating sign matrices is given. The derivation is based on the recent results of 
Stroganov for the corresponding generating function. As a result, A(n, r;3)'s are rep- 
resented as 1-fold sums which can also be written in terms of terminating 4F3 series of 
argument 1/4. 



1. Introduction 

An alternating sign matrix (ASM) is a matrix of l's, O's and — l's such that in each row and in 
each column the first and the last nonzero entry is 1, and all nonzero entries alternate in sign. 
There are many nice results concerning ASMs (for a review, see [1]). The most celebrated 
result gives the total number A(n) ofnxn ASMs. It was conjectured by Robbins, Mills and 
Rumsey [2, 3] and subsequently proved by Zeilberger [4] and Kuperberg [5] that 

... A (3*-2)!(fc-l)! n (3k - 2)\ 

^»)=n (a _ 1)l(2i _ 2)l =n(srr^- o-« 

Yet another formula was conjectured in [2, 3], and proved in [6], concerning the refined enu- 
meration of ASMs, namely, that the number of n x n ASMs with their sole 1 of the first (or 
last) column (or row) at the r-th entry, denoted as A(n, r), is 



\ n- 1 I \ n- 1 ) 
/3n-2\ 
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A(n,r)= K " r 1 A (n). (1.2) 



One more conjecture of Robbins, Mills and Rumsey, proved by Kuperberg [5], gives the num- 
ber of 3-enumerated ASMs. In general, ^-enumeration counts ASMs with a weight x k , where 
k is the total number of —l's in each matrix (the number x here should not be confused with 
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the variable x below). The total number of 3-enumerated ASMs is denoted as A(n; 3) and the 
result reads 



A(2m+1;3) = 3" ,( "" n ^[ 
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,(3m + 2)\m\ 



A(2m+2; 3) = 3"^ '— ¥ A(2m+1; 3) . 

,^L(m + *)!j [(2m + I)!] 2 

(1.3) 

In the present paper we study the problem of the refined 3-enumeration of ASMs, i.e., we 
derive the numbers A(n, r; 3). The main result can be summarized as follows. 

Theorem (Main result). The refined 3-enumeration of ASMs, A(n, r; 3), is given by 



b(m, r — 1) + b(m, r — 2) 
A(2m + 2, r; 3) = • ' V ' - A(2m + 2; 3) 

2 6(m,r-l) + 56(m,r-2) + 26(m,r-3) 

v4(2m + 3, r; 3) = A(2m + 3; 3) 
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where the quantities b(m, a), obeying b(m, a) = b(m, 2m — a), are given by 

[a/2] 



(1.4) 
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i=max(0,a-m) 
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/or a = 0, 1, ... , 2m, while they are assumed to be zero for all other values of a. 

As a comment to the main result it is to be mentioned that the numbers A(n, r; 3) appear not 
to be writable as a single hypergeometric term, i.e., not to be 'round' (or 'smooth'), contrarily 
to the expressions for A(n), A(n, r) and A(n; 3). For instance, direct inspection of the results 
of computer experiments for A(n,r; 3) shows the appearance of large prime factors in their 
prime factorizations, strongly suggesting that no answer for A(n, r; 3) can be given in the form 
of a single hypergeometric term. This explains why no conjecture concerning A(n, r; 3) was 
given previously, and it might moreover imply that any closed expression for A(n, r; 3) could 
be at best a sum of several 'round' terms. In view of this, the expression in formulae (11.41) and 
(11.51) . even if not as elegant and neat as other results in the context of ASM enumerations, is 
probably the best that can be achieved: a 1-fold sum of 'round' terms. 

Our proof of the Theorem is based on a certain representation for the generating function 
for the numbers A(n, r; 3), recently obtained by Stroganov [7]. In the next Section we recall 
the relevant formulae of paper [7] and shortly review their backgrounds. In Section 3 we apply 
standard relations for Gauss hypergeometric function to simplify significantly Stroganov's ex- 
pression for the generating function of 3-enumerated ASMs. In Section 4, from this simplified 
expression we obtain explicit formulae for A(n, r; 3), thus proving the Theorem. 
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2. Preliminaries 



In studying the refined 3-enumeration ofnxn ASMs it is convenient to define the generating 
function 



I J^^A(n,r;3)r- 1 , lf< 3 >(l) = 1 . (2.1) 

A simple property of the generating function is Hn\t) = t 71 " 1 H^^' 1 ) , which follows from 
the relation A{n, r; 3) = A(n, n — r + 1; 3) expressing left-right (top-bottom) symmetry within 
the set ofnxn ASMs. 

In paper [7] the following defining properties of the generating function have been estab- 
lished. The generating functions Hn\t) for n even, n = 2m + 2, and for n odd, n = 2m + 3, 
are related by the formulae 

flSU) = ^ ft»W • "£«M = (2< + y + 2> . (2-2) 

where B 2m (t) is some polynomial of degree 2m. It is convenient to define E m (t) := t~ m B 2m (t) 
so that i? m (t) = E m (t~ l ). The function ^(t) can be related to another function V m (x), 
obeying, in turn, the relation V m (x) = V rn (x~ l ), by means of the transformation 

E m(t) = . V ;i X) 1]m , t = (2.3) 

(x — 1 + x L ) m qx — 1 

or, inversely, 

3 m E m (t) t + g 

Knlx) = ; -— , X = - (2.4) 

v ; (t + l + t- l ) m ' tg+ 1 

where (and everywhere below) the following shorthand notation is used 

q := exp(i7r/3) . (2.5) 

The function 

h m (x) := (x - x- l ) 2m+ \x + 2 + x- 1 ) V m (x) (2.6) 

is found to be 

(3m + 2 \ 

9m(x) + - f m (x) j (2.7) 

where the functions f m (x) and g m (x) are given by 

g m (x) ■= £ + fc 2/3 ) f ) f^ 2-6 * - 1- 3 — 2+6 ') (2.8) 

/-(*) := E + ,. 1/3 ) (" "-f ) ( l3 ™ +1 "' - I "™" +6 ') C ' 9) 

fc=0 ^ ' ^ ' 
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and c m is some constant such that V m (l) = E m (l) = 1. For the reader's convenience we note 
that the variable x here and the variable u of Ref. [7] are related by e 2m = — x^ 1 . 

The result just described for the generating function of 3-enumerated ASMs was obtained 
by analyzing Izergin's determinant formula [8] for the partition function of the six-vertex model 
with domain wall boundary conditions (also known as 'square ice'). The one-to-one correspon- 
dence of square ice states with ASMs was pointed out in [9] (see also [10]) and fruitfully applied 
by Kuperberg in [5]. The approach proposed in [11] for 1 -enumerated ASMs and extended to 
3-enumerated ASMs in [7], uses a certain functional equation for the square ice partition func- 
tion. This functional equation is nothing but the Baxter's T-Q equation for the ground state of 
XXZ Heisenberg spin- 1/2 chain at A = —1/2 and with odd number of sites iV = 2m + 1, 
studied previously in [12]. 

In particular, expression (12.71) for h m (x) arises as the solution of the functional equation 



which is supplemented by the condition that its solution y(x) = h m (x) must have the structure 
(12.61) where the function V m (x) is supposed to be such that V m (e lip ) is an even trigonometric 
polynomial of degree m in the variable if. This condition allows one to reduce the functional 
equation to a finite set of linear algebraic equations. These equations can be solved explicitly 
thus leading to the formula for h m (x). 

The functions f m (x) and g m (x) are also solutions of equation (12. 1 Ob . and by construction 
they have the structure 



f m (x) = (x-x- 1 ) 2m+1 g m (x) , g m (x) = {x-x- x f m+1 P m+l {x), (2.11) 



where Q m (x) and P m +i (x) are such that Q m (e lv ) and P m +\ (e lip ) are even trigonometric polyno- 
mials of degrees m and m + 1 respectively. These two functions are two independent solutions 
of T-Q equation. Note that as the second independent solution (i.e., that solution which is of 



of the form P m+ i(x) + a Q m (x), see [13]. Thus, (12.71) just implies that 3-enumerated ASMs 
are described by the second independent solution of T-Q equation. The first solution, Q m (x), 
was shown in [1 1] to be related to 1 -enumerated ASMs. 

Namely, in the case of 1 -enumerated ASMs it was shown that 



y(x) + y(q 2 x) + y(q 2 x) = 



(2.10) 



higher degree) one can also regard the function P m +\ {%) + 



3m+2 
3m+l 



Q m (x) as well as any function 





qx — q L x 



— q x 



(2.12) 



where the generating function Hn if) is defined similarly to (12.11) by 




r=l 




(2.13) 
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To reproduce formulae (ll.lt and (11.21) from (12.121) it was proposed to use the differential equa- 
tion satisfied by f m (x), which reads 



[dl - 6m cot 3^ d v - (3m + l)(3m - 1)] / m (e il ^) = 



(2.14) 



Fortunately, the corresponding equation on the generating function appears to be a hypergeo- 
metric equation 



[{t-l)td 2 t +2(t + n-l)d t -n{n-l)]H^\t) = 

and its polynomial solution gives 

(2n- l)!(2n-2) 



iFi 



—n + 1 , n 
-2n + 2 



(3n - 2)! (n - 1)! 

Here the proper normalization, see (I2.13I) . follows from the Chu-Vandermonde identity 



(2.15) 



(2.16) 



—m , b 

c 



{c-b) t 

(c)m 



(a) m := a(a + 1) • • • (a + m — 1). 



(2.17) 



Since A(n, 1) = A(n — 1), Eqn. (12.161) results in an elementary recursion 



A(n 



(2n - 1)! (2n - 2)! 



(2.18) 



A(n) n v ' (3n-2)!(n-l)! 

which, supplemented by the initial condition A(l) = 1, gives formula (ll.lt . Expanding RHS 
of (I2.16t in power series in t and multiplying these coefficients by A(n), one reproduces, in 
turn, formula (II .21) . 

Unfortunately, a similar procedure cannot be applied in the case of 3-enumerated ASMs. It 
can be verified directly that the function h m (x) solves the equation 



si 



(2m + 1) cot 3<p 



sin3<£> 



- (3m + l)(3m + 2) U m (e i¥3 ) = 0. (2.19) 



This equation leads to a rather bulky non-hypergeometric equation for the generating function 
H { n\t), which can hardly be solved by standard means. Moreover, even though the function 
Q m (x) can easily be deduced from the equation for f m (x) above, the differential equation for 
9m(x) (for this equation, see [12]) does not allow one to find the function P m+ i(x) in a similar 
way. As a matter of fact, only an implicit answer for Hn(t) in terms of two recurrences, gen- 
erated by three-terms recurrences for f m (x) and g m (x), was given in paper [7]. Here we would 
like just to mention that from differential equation (I2.19t it can easily be seen that the opera- 
tor (sin 3(p)~ 1 d lfi plays the role of a 'lowering' operator in the set of trigonometric polynomials 
{h m (e 1,p )}. It is thus straightforward to rewrite Eqn. ( I2.19t as a three term recurrence for h m (x) 
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and, hence, to obtain a single recurrence for Hn (t). However this point of view will not be 
developed further, since it amounts merely to a reformulation of the problem. Instead, in the 
next Section we show how to overcome the difficulty and find the function Hn^ (t) explicitly. 

As a last comment here we would like to mention that, although the previous results on the 
refined 3-enumeration of ASMs run over papers [7, 11, 12], the present research have originated 
from our study of the square ice boundary correlators obtained in [14]. Within our study the 
refined enumerations of ASMs arise as solutions of some differential equations rather than 
functional ones, in particular, the expression for h m (x) arises as a solution of (I2.19I ). The 
illustration of this approach is out of the scope of the present paper and will be given elsewhere 
[15]; here we restrict ourselves to obtaining the numbers A(n, r; 3) from the known h m (x). 



3. The generating function 

We begin with noticing that the functions f m (x) and g m (x) can be also written as follows 



(X 



T(m+ 1/3) 
m!r(l/3) 



r(m + 4/3) 
m!r(4/3) 



x 



3m+2 771 
2^1 



-m , — m — 2/3 
1/3 



x 



X 



-3m-2 



x 



-3m+l 771 
2-^1 



— m , — m — 2/3 
1/3 

— m , — m + 1/3 
4/3 

— m , — m + 1/3 
4/3 



x 



(3.1) 



x 



— X 



3m— 1 



x 



-6 



(3.2) 



Since the parameters of the hypergeometric functions entering these expressions differ by in- 
tegers one can expect that g m (x) and f m (x) are connected by some three-term relations via 
Gauss relations (see, e.g., Sect. 2.8 of [16]). Indeed, using Gauss relations it can be shown that 



2^1 



-m , — m — 2/3 
1/3 



3m + 4 



2*1 



m — 1 , — m — 2/3 

4/3 

3m + 2 



—m , — m + 1/3 
4/3 



z) (3.3) 



and therefore we have the following relation 

3m + 2 



9m\X) 



2(3m+ 1) 
Substituting (l3~4b into (|277l) we obtain 

3m + 2 



(x 3 + x 3 ) f m (x) - | f™ + 1 \ f m +i{x). 



2(3m+ 1) 



(3.4) 



2(3m+ 1) 



3 777, — (— 3 

(X 3 + 2 + X~ 3 ) f rn (x) - — f m +l(x) 

3m + 2 



(3.5) 
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Comparing Eqns. (13.41) and (13.51) with Eqns. (12.61) and (12.111) we find that 



*m+l(+ 

and 



3m + 2 
~ 2(3m+ 1) 

3m + 2 



(x 3 + x 3 )Q m (x) 



2(3m+ 1) 



X — 1 + X 



-1\2 



3(m + 1) 
2(3m+ 1) 

3m + 3 , 



(x-x 1 ) 2 Q m+ i(x) 



X-2 + X ^Qm+iW 



(3.6) 



(3.7) 



3m + 2 

Hence all functions in question are expressed in terms of Q m (xYs. 

The advantage of this approach is based on the fact that for the function Q m (x) the follow 
ing explicit formula can be found 



(2m) 



[ml 



qx~ 



q 1 x 



q-q 1 



2*1 



— m , m + 1 
—2m 



qx — q 1 x 1 



qx 1 — q 1 x 



(3.8) 



As already mentioned in our preliminary comments, the function Q m (x) can be deduced from 
the differential equation for f m (x) above (the proper normalization, for instance, can be found 
using the three-term relation for / m (x)'s). Here we give a proof of (13.81) by a straightforward 



transformation of the function f m (x) to the form f m (x) 



(x — x 



-l\2m+l 



The key identity which is to be used here is the so-called cubic transformation of Gauss 
hypergeometric function [16] which in its most symmetric form reads: 



r(q) 
r(2/3) 



2F1 



a + 1/3 , a 
2/3 

_ g-3a+l 



T(a + 2/3) 

u z r(4/3) 2jPi 

I - z\-3« T(3a) 



1-a; 



2-^1 



a + 1/3, a + 2/3 
4/3 
a + 1/3 , 3a 
2a + 2/3 



z — C<J 



UJ- 



I- (3.9) 



T(2a + 2/3) 

Here a; is a primitive cubic root of unity, uj = exp(±2i7r/3), and a is arbitrary parameter. To 
show that indeed the cubic transformation is relevant to our case, let us rewrite (13.21) in the form 
consistent with LHS of (13.91) . Taking into account that 



2*1 



— m , — m + 1/3 
4/3 



r(i/3)r(4/3) 



T(-m + l/3)r(m + 4/3) 



-z) 2 



Ft 



-m , — m — 1/3 
2/3 



and 



r(i/3) 

T(m + 4/3) 



-1 



>m+l 



r(-m - 1/3) 



T(m + 4/3) 



„(3m + r 



T(2/3) ' T(-m + 1/3) 
it is easy to see that (13.21) can be rewritten in the form 

^-l) m+1 (3m + 1 



3 3m+l m \ 



(3.10) 
(3.11) 



fm(x) 



33m+l 



X 



,3m+l 



ml 



T(-m- 1/3) 

. r(2/3) 



2*1 



r(-m + l/3) 

+ 3^ — , — 2^1 



r(4/3) 



— m , — m — 1/3 
2/3 

— m , — m + 1/3 
4/3 



x 



(3.12) 
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Clearly, both terms in the brackets are the same as in LHS of ( 13.91) provided the parameter a is 
specialized to the value a = — m — 1/3. 

To apply the cubic transformation to (13.121) we first define 



W{a;z) :-- 



so that 



r(q) 
r(2/3) 



2^1 



a + 1/3 , a 
2/3 



z 6 + z 



r(g + 2/3) 
r(4/3) 



2^1 



a + 1/3, a + 2/3 
4/3 



fm(x) 



-1 



im+1 



(3m + 1)! 



x 



3m+l 



3 3m+l ( m l)2 

Next, we note that for a sum of two terms one can always write 

q 



W(-m- 1/3; x- 2 ). 



X + Y = —^— 1 (X-q- 2 Y) 



Q 



(3.13) 
(3.14) 

(3.15) 



g — q x ' q — q~ 

and if g = exp(i7r/3), which is exactly the case, one can set to = q 2 for the first pair of terms 
and uj = q~ 2 for the second one. This recipe allows one to apply the cubic transformation, that 
gives 



W(a;z) 



)-3o+l 



r(3o) 



:i-z) 



-3a 



T(2a + 2/3) q(l - g2)-3a+i _ 



2^1 



a + 1/3 , 3a 
2a + 2/3 



q 



l-z 



+ q 



3a+l 



•ft 



a + 1/3 , 3a 
2a + 2/3 



-2 Z ~ (i 



1 - z 



To obtain a new formula for f m (x) via (13.141) we have to evaluate now the limit a 
of (I3.16I) . The limit of the pre-factor can be easily found due to 

lim T(y3a ^ 2 ( _1 ) m+1 ( 2m ) ! 

a _>_™ 1/3 r(2a + 2/3) ~ 3 (3m + 1)! ' 

To find the limit of the expression in the brackets in (13.161) we note that 



2 Z ~T 

q =: s 

1 — z 

and hence the following formula can be used 



q 



-2 z - q~ 



l-z 



1-8 



. (3.16) 

-m — 1/3 

(3.17) 
(3.18) 



lim 



,'a + l, 3a 
' 2a + f 



s \+q 



-.',„-;- ! / ■ / a + 3 ' 3a 



1-S 



2^1 



-m , —3m — 1 
—2m 

(3.19) 

Formula ( 13.191 ) can be proved, for instance, by virtue of standard analytic continuation formulae 
for the hypergeometric function (see, e.g., Eqns. (1) and (2) in §2.10 of [16]). Collecting 
formulae we arrive to the expression 



W{-m- 1/3; z) 



j2m+l 



(2m)! (l-z 



,3m+l 



(3m +1)! (q-q- 1 )' 



—m , —3m — 1 
—2m 



1 - z 



q 



. (3.20) 
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Finally, substituting this expression into (13.141) and using the identity 

' a , c — b z 



a , b 
c 



(3.21) 



we obtain 



fm{x) 



(2m) 



(x — X 



—1 _i \ m 



[m\ 



q-q- 1 



Fx 



—m , m + 1 
—2m 



qx — q l x 1 



gar 



(3.22) 



Obviously, this expression leads directly to (13.81) which is thus proved. 

Hence, we have just shown that a particular solution of equation (12. 1 Ob . the function f m (x), 
Eqn. (12.9b . is connected to the 'first' solution of the Baxter T-Q equation, Q m (x), Eqn. (I3.8I) . 
via the cubic transformation. Formulae (13.41) and (13.51) mean that the same transformation is 
responsible for relationship of g m (x) and h m (x) with the second independent solution of T- 
Q equation, which can be chosen either as P m+ i(x), in the case of g m (x), or as P m+ x(x) + 
§23 QmO) = + 2 + x' 1 ) V m {x), in the case of h m (x). 

To write down the resulting expressions for functions P m+ i(x) and V m (x), and for the 
purpose of subsequent transformation of these expressions, we define 



$( fc )i 



x) 



—l — i \ m 
qx — q £ 



q-q 



-i 



Ft 



— m , k + 1 
— m — A; 



qx — q l x 1 



1 — g 1 x 



(3.23) 



Note that for positive integer m, function $^ (x) possesses the property $^ (x) 
and thus, obviously, $^ (x) is a polynomial of degree m in the variable u := x + x~ x . 
In particular, Eqn. (13.81) now reads 

(2m)! 



Qm [X 

Formulae (13.61) and (13.71) give 

3m + 2 (2m)! 



3 m (m!) 



2 * m 



X) 



(3.24) 



•PmH-l(z) 



2(3m+ l)3 m (m!) 2 _ 



(x 3 + 2T 3 )$L m) ( 



2(2W+l) ^-,- 1 ) 2 e + 1 1) (-) 



(2m)! (2m + 1)! 
m! (3m + 1)! 



3m + 2 
2(2m+ 1 



(3.25) 



-1\ ifXm+X). 



3m + 2 ^-2 + 0^ + TW 



(3.26) 



Here the expression for V m (x) is written according to the proper normalization of this function, 
V m (l) = 1. The normalization can be verified by virtue of Chu-Vandermonde identity (12.171) 
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and it corresponds to the choice 



c m = (3m + 1] 



3 m+1 m! (2m + 2)! 



(3.27) 



(3m + 3)! 

Formulae (13.251) and (13.261) for functions P m+1 (x) and V m (x) are not the final expressions 
for them yet since they can be notably simplified. Indeed, again using Gauss relations one can 
prove, for instance, the following relations for the functions (13.231) 



, \ . i \ fh\ / \ m(m + 2k + 1) . 9 

= iX + * > *» <*> - 3(m + fc + l)(m + fe ) ^ + 1 + 



(3.28) 
(3.29) 



2(m + fc + 1) ~ m v ~ y ' 2(m + k + 1 

These two relations can be regarded as basic relations; many other tree-term relations connect- 
ing $m (x)'s can be derived from (13.281) and (13.291) . Using these three-term relations we find 
that for the function P m +i(x), in particular, the following formula is valid 



Pm+l \ x ) 



(2m)! 



3 m m\(m + 1 
Similarly, for the function V m (x) we obtain 
(2m)! (2m + 2)! 



(3m + 2) ^Z^\x) - (2m + 1) $ 



(m) , 
m+l' 



(m + 1)! (3m + 2)! 



(2m + 1) -m{x-l + aT 1 ) $ 



m—l 



(3.30) 



(3.31) 



This formula for V m (x) can be further simplified, due to (13.291) . to be written solely in terms of 
$ m (/c)'s, but (13.311) is more convenient for the purpose of obtaining the function E m {t), which 
is of primary interest for 3-enumerated ASMs. Applying (12.41) we obtain 



E m (t) 



(2m)! (2m + 2)! 
3 m (m + 1)! (3m + 2)! 



(2m + l)(t + 2) m 2 F : 



-m , m + 2 



- 3m (t + 2) m - x 2 Fi 



-2m - 1 
m + 1 , m + 2 
—2m 



2t 



t{t + 2) 
1 + 2* ' 



t(t + 2) 

Formulae (13.30b . (13.311) and, especially, (13.321) are the main results of this Section. 



(3.32) 



4. The numbers A(n, r; 3) 

Let us begin with showing how Eqn. (11.31) for A(n; 3) is recovered from the just obtained 
expression for the generating function. Here the property A(n, 1; 3) = A(n — 1; 3) is to be 
used. It implies 



A(2m + 2;3) 2m+2V ; 2 V ; ' A(2m + 3; 3) 2m+3V ; 9 V ; 

(4.1) 
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Hence, one has the recurrences 

A(2m + 3;3) = =A(2m+l:3), A(l:3) = l; (4.2) 

[B 2m (o)] 2 

A(2m + 2; 3) = j? ^A(2m; 3), A(2; 3) = 2. (4.3) 

Recall that B 2m {t) = t m E m {t) in (l2~2l . From (l3~32l one finds 

(2m + 1)! (2m + 2)! 
B2m(0) " 3-(m + l)!(3m + 2)!- (44) 

Substituting (14.41) into (14.21) and (14.31) . and solving the recurrences one obtains (11.31) . 

Let us now turn to our main target, A(n, r; 3), which describes the refined 3-enumeration 
of ASMs. Formulae (12.21) imply that 



b(m, r — 1) + b(m, r — 2 ) 
A{2m + 2, r; 3) = ^ V ' - A(2m + 2; 3) 

(4.5) 

2 b(m, r — 1) + 5 b(m, r — 2) + 2 6(m, r — 3) 

A{2m + 3, r; 3) = — K -^- — — A(2m + 3; 3) 

9 

where b(m, a) are defined as 



2 m 



B 2m (t) = t m E m (t) =: Km, a) t a (4.6) 



Q = 



and are assumed to be zero if a is out of the range of values 0, 1, ... , 2m. To find these 
coefficients we expand (13.321) in power series in t, thus expressing B 2m (t) as a triple sum. 
Chu- Vandermonde summation formula can now be applied to the sum with respect to the index 
defining the hypergeometric series in (13.32b . thus expressing B 2m (t) as a double sum. These 
two summations can be rearranged in such a way that one of them becomes with respect to a 
while the other one defines the coefficients of power expansion in t. We obtain 



. 2m + l !m! [ M . n . /3m + 3\ 

b(m,a ) = —, ^7 > 2m + 2-a + 2f) 

v ' ; 3 m 3m + 2 ! ^ v \a-2i) 

' «=max(0,Q-m) 



x 2m + l-a + l\ m + l + l\ 2e (4?) 
\ m + 1 J \ m + 1 J 

Here [a/2] denotes integer part of a/2. Formulae (14.51) and (14.71) complete the proof of the 
Theorem, which constitutes the main result of the present paper. Here we would like to men- 
tion that in terms of terminating hypergeometric series the last formula may be rewritten, for 
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instance, as follows 



2«('3m+3\ /2m+l-a\ 
7 I \ V a I V m+1 / 

m, « = -77, -777 

6 \ m+1 ) 



2 4^3 



-(a- l)/2, -a/2, m + 2, 2m + 2 - a 
(3m + 4 - a)/2 , (3m + 5 - a)/2 , m - a + 1 



a , / -(a- l)/2, -a/2 + 1, m + 2, 2m + 2 - a 
m+1 4 3 V(3m + 4-a)/2, (3m + 5-a)/2, m-a + 1 



y 

(4.8) 



This formula is valid for a — 0, 1, . . . , m (a similar expression for a = m + 1, m + 2, . . . , 2m 
can be simply obtained through the replacement a — ► 2m — a in RHS of (14.81) ). The two 4 F 3 
in (14.81) can be further combined into a single 5F4 . Analogous formulae for b(m, a) in terms 
of terminating hypergeometric series of argument 4 may be written down as well. Analyzing 
these expressions we were anyway unable to perform the sum in (14.71) in a closed form, even if 
very suggestive similarities can be found with known summation formulae, see §§7.5 and 7.6, 
especially §7.6.4, of Ref. [17]. We therefore regard Eqns. (14.51) and (14.71) as the most compact 
formulae for the numbers A(n,r;3) available at the present moment. 

To conclude, let us comment the just obtained result in application to the square ice. The 
quantity of interest here is nA(n, r; 3)/A(n; 3) which plays the role of spatial derivative of 
the boundary polarization of the square ice with a particular choice of the vertex weights. It 
is interesting to study the behavior of this quantity in the large n and the large r limits, with 
the ratio £ = r/n kept fixed, < £ < 1. In the context of the square ice this scaling limit 
corresponds to the continuous limit, with \ jn playing the role of a lattice spacing. For the 
square ice in the disordered regime (which corresponds to x-enumerations with < x < 4) a 
Heaviside step-function behavior #(£ — 1/2) of the boundary polarization is expected [14]. In 
ASM enumeration language this corresponds to 

lim n A y; r]X ) =8{£-l/2) for < x < 4. (4.9) 
n,r^oo A(n:x) v ' ' 

r/n=5 

For instance, in the case of 1-enumeration of ASMs (i.e., for x — 1) the validity of this formula 
can be verified directly from (11.21) by virtue of Stirling formula. This result is also valid for 
the simple case of 2-enumerated ASMs whose refined enumeration is just A(n, r; 2) /A(n; 2) = 
/2 n_1 ; for the result for the whole free-fermion line see [14]. Formulae (14.51) and (14.71) 
allows one to study the scaling limit of the expression nA(n, r; 3)/A(n; 3) as well. In this limit 
the sum in (14.71) turns into an integral and standard saddle point method can be applied. In this 
way we find 

lim n ^ ,r;3) 1/2). (4.10) 

r/n=5 

This confirms (14.91) also in the x = 3 case. 
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